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ABSTRACT 



We study gravitino dark matter and slow gravitino decays within the framework of 
R-violating supersymmetry, with particular emphasis on the flavour dependence of the 
branching ratios and the allowed R-violating couplings. The dominant decay modes and 
final state products turn out to be very sensitive to the R-violating hierarchies. Mixing ef- 
fects can be crucial in correctly deriving the relative magnitude of the various contributions, 
particularly for heavy flavours with phase space suppression. The study of the strength of 
different decay rates for the gravitino is also correlated to collider signatures expected from 
decays of the Next-to-Lightest Supersymmetric Particle (NLSP) and to single superparticle 
production. 
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1 Introduction 



Recently, there has been renewed interest in the possibihty of having gravitino dark matter 
within the framework of R- violating supersymmetry [1,2], which occurs if the gravitino 
decays are slow enough for its lifetime to be larger than the age of the universe [3,4]. This 
is an exciting possibility that allows supersymmetric dark matter, even if the symmetries of 
the fundamental theory result in an unstable Lightest Supersymmetric Particle (LSP) [5-7]. 

This is what happens if, in addition to the couplings that generate the fermion and 
Higgs masses 

(XH1H2 + m^UEjHi + m'^QiDjHi + m''QiEjH2, (1.1) 
we also have R-violating couplings of the form 

hLiH2 + XLiLjEk + X'LiQjDk + X"UiD,Dk. (1.2) 

In the above, Hi^2 are the Higgs superfields, L{Q) are the left-handed lepton (quark) doublet 
superfields, and E {D, U) are the corresponding left-handed singlet fields. The first three 
couplings in (11. 2p violate lepton number, while the fourth violates baryon number. 

The stricter bounds on R-violating operators come from proton stability, and R-parity 
[8], which forbids all lepton and baryon number violating operators, is one of the possible 
solutions. However, this is not the only symmetry that can guarantee proton stability; 
baryon or lepton parities [9, 10] can also exclude the simultaneous presence of dangerous 
LQD and UDD couplings [11]. Experimental constraints from the non-observation of 
modifications to Standard Model rates, or of possible exotic processes [12] also impose 
additional bound^. 

R-violating supersymmetry results in a very rich phenomenology. In the presence of the 
additional operators, the NLSP can decay into conventional Standard Model particles. The 
missing energy signature of the Minimal Supersymmetric Standard Model (MSSM) [14] is 
substituted by multi-lepton and/or multi-jet events. In addition to the consequences for 
collider searches, R- violation implies that any gravitinos that have been thermally produced 
after a period of inflation, are also unstable. 

Gravitinos have three main decay modes: via tree-level three-body decays to fermions [4] , 
via two-body decays to neutrino and photon due to neutrino-neutralino mixing [1,3], and 
via one-loop decays to neutrino and photon, generated by the trilinear couplings [2] q In 
all three cases, the very large suppression 1/Mp of the gravitino vertex, where Mp is the 
reduced Planck scale, plus additional suppression from phase space, mixing and loop fac- 
tors, respectively, result in large gravitino lifetimes. For a wide set of R-parity violating 
couplings and gravitino masses these exceed the age of the universe. Moreover, the photon 
flux from these decays could be able to explain the apparent excess in the extragalactic 
diffuse gamma-ray flux in the re-analysis of the EGRET data [15, 16]. 

^Additional strong constraints can be obtained from the observation of NLSP decays to a gravitino LSP, 
with a photon or lepton plus missing energy signature [13]. 

^For heavy gravitinos, there is also the possibility of producing massive gauge bosons. However, for 
trilinear couplings and the range of parameters considered here, these contributions are subdominant. 
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The branching ratios for gravitino decays are sensitive to the flavour structure of the 
R-violating operators. In the case of O (GeV) gravitinos, the presence of tau or bottom 
quarks in the final state significantly enhances the branching ratio of radiative decays with 
respect to the tree-level ones, while for "super-light" gravitinos, as in [17], gravitinos are 
essentially stable with respect to the three-body decays. Moreover, in the case of non-zero 
\"lJ^DjDk only — with a top quark final state — gravitinos lighter than rrit have a maximal 
stabihty, modulo mixing effects, which we will discuss in a subsequent sectioij^l. 

In [2], gravitino decays were studied for LLE operators that give rise to both loop and 
tree-level decays, with a tau or a muon in the loop. Here, we extend the results in the 
following way: 

(i) We look at flavour effects in more detail, making the link with fermion mass hier- 
archies. Within this framework we comment on the relative magnitudes for bilinear 
and trilinear R-violation and what are the implications for gravitino decays. 

(ii) We extend the discussion to all 45 LLE, LQD and UDD operators, paying particular 
attention to the different features of the various decay modes and possible bounds 
from gamma-ray measurements. 

(iii) We consider possible implications of mixing effects, which in certain cases can be 
quite significant. For instance, for the UsDjDk operator, the expected decay depends 
very sensitively on the right quark mixing (for which little information is available). 

(iv) We link the above with probes of R-parity violation at the LHC, in particularly NLSP 
decays, which may yield interesting signatures. 

We begin in Section [2] by describing the various modes of gravitino decays with trilinear 
couplings and the calculation of the resulting extragalactic diffuse photon flux. In Section [3] 
we discuss possible flavour structures for R-parity violating operators, before we look at 
the consequences for gravitino decays in Section HI with particular attention to bounds 
from gamma-ray measurements. We continue with the corresponding prospects for hadron 
colliders in Section [5], before concluding in Section [6l 

2 Gravitino Decays 

As already discussed, trihnear R-violating operators may cause gravitinos to decay via two 
different channels: 

• Via two-body radiative loop decays to neutrino and photon (Fig. [1]) [2]. 

• Via tree- level decays to fermions (Fig. [2]) [4] . 



•^For an operator of the form X'LiQ^Dk this argument does not hold, since, when we pass from superfields 
to component fields the LiQ^ part can become £it or Uib. 
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Figure 1: Basic set of Feynman diagrams for radiative gravitino decay, shown for (s)fermion 
loops. In the case of (s)quarks, the neutrino is coupled to down-type quark to preserve SU(2) 
invariance. Arrows denote flow of fermion number for left- chiral fields. 




Figure 2: Three-body decay of a gravitino via an R-parity violating coupling. There are 
three contributing diagrams where the sfermion carries any one of the three indices i, j and 
k of the corresponding operators. 



The decay rates have been presented in detail in the original references, and for com- 
pleteness are briefly summarised in Appendices A and B, respectively. For light gravitino 
masses and appropriate fermions in the loop the radiative decays may dominate. Indeed, 
as we shall see, even when the three-body decay involving an intermediate sfermion / is 
well above the kinematical threshold at 2mj, the radiative dominance is still present. The 
behaviour of the decay rates is controlled by the mass dependence of the decay width: for 
the three-body decay Tq oc m^, while for the radiative decay F^^ oc at low gravitino 
masses. The latter occurs since the gravitational coupling compensates for the relatively 
high loop mass by its increasing strength for higher loop momenta. Because of the helicity 
structure of the couplings, the two-body decay width is also strongly dependent on the mass 
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of the fermion in the loop, oc mf at low gravitino masses, implying significantly shorter 
lifetimes for dominant third generation couplings. 

To constitute a realistic dark matter candidate, the gravitino lifetime should exceed the 
age of the universe. Moreover, the photon flux from gravitino decays has to be consistent 
with observations. The diffuse extra-galactic gamma ray flux of energy E from the gravitino 
decays is described by a integral over red-shift z given by [18] 

^(^).^.i^.2^c„r..-^^f:L=, ,2.1) 



where y = z and dN^/dE is the gamma ray spectrum from the gravitino decay. Here 

C, = ^^P' and K = ^. (2.2) 

For the radiative gravitino decay dN^/dE = 6{E — mQ/2) and Eq. (12. ip simplifies to [1] 

F{E) = E^^ = BR{G ^ 7z/)C^(l + Kx^y^/^x''/^e{l - x), (2.3) 
dE 

where x = 2E/mQ. In the case of three-body decays the hadronization of the produced 
particles and the resulting photon spectrum have been calculated using PYTHIA 6.4 [19]. 
The photons from the three-body decays come mostly from internal bremsstrahlung off 
leptons and from 7r° decays. 

Using the original EGRET analysis [15], with a power law description of the extragalac- 
tic flux as 

= 1.37 • 10"^ ( ^-^) GeV cm'^ sr-^ s^', (2.4) 
dE \ E J 

in the energy range 30 MeV to 100 GeV, severe bounds on gravitino decays and thus on 
the allowed combinations of gravitino masses and R-violating couplings can be derived. 
For comparison, predictions for photonic spectra from gravitino decays through neutrino- 
neutralino mixing, and also possible antimatter signatures of gravitino dark matter, have 
recently been studied in [18] and [20]. 



3 Flavour Structure and Hierarchies of R-violating 
Operators 

The implication of radiative gravitino decays as compared to the tree-level ones, clearly de- 
pends on the flavour structure of the R-violating operators involved. For higher generations 
the radiative decay widths become larger and the tree-level diagrams suppressed due to lim- 
ited phase space. Most phenomenological studies assume a single operator-dominance. This 
can be motivated by the fact that the Yukawa couplings that generate fermion masses also 
have large hierarchies. However, in principle, one may try to relate R-violating hierarchies 
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to those of fermion masses [21,22], using models with family symmetries. When exact, 
the latter allow only the third generation fermions to become massive, while the remaining 
masses are generated by the spontaneous breaking of this symmetry (see below). If R parity 
is violated, couplings with different family charges will also appear with different powers of 
the family symmetry-breaking parameter, and thus with different magnitudes. 

Moreover, one would have to appropriately take into account mixing effects. Indeed, 
even with the common assumption of single /^-violating operator dominance, this would 
be true only for the basis of current eigenstates for quarks and leptons, while, in the mass- 
eigenstate basis, there would be several operators corresponding to the original dominant 
one in the current basis. In addition, the fact that there are strict bounds on some operators, 
implies that mixing effects may in given models generate additional bounds on couplings 
that at a first glance look less constrained. This has been analysed in detail in [22] , where 
it was shown that in theories with strong correlations between operators (such as left-right 
symmetric models), the effects can be particularly significant. 

The starting point in such considerations, is to assume a f/(l) flavour symmetry, with 
the charges of the Standard Model fields denoted as in Table [TJ 





Qi Ui D, Li Ei N H2 Hi 


U(l) 


ai f3i 7i Ci di d -03 - /^s w 



Table 1: Notation for possible U(l) charges of the various Standard Model fields, where i 
is a generation index. 

The flavour charge of H2 is chosen so that the operator that generates the top quark 
mass {Q3U3H2) has a zero U(l) flavour charge and thus is allowed at zeroth order, as it 
should be, since the top quark is significantly heavier than the rest. The remaining matrix 
elements may be generated when the U(l) symmetry is spontaneously broken [23,24] by 
fields 9, 9 that are singlets of the Standard Model gauge group, with U(l) charges that are 
in most cases taken to be ±1, respectively. For instance, for a^ = f3i and |a3 — = ±1 as 
in [9], the charm mass comes about by a term Q2U2H2 {{9)/M) or Q2U2H2 {{9)/M) where 
M is the heavy scale of the theory. 

One may generalise the above to non-abelian flavour symmetries, and, as an example, 
the following mass matrices have been proposed [25] : 

/ e3 \ / e3 / \ 

M^P oc e3 e2 e2 , M'^^^'^ oc eM , oc . 

V e2 1 / V e3 1 / V 1 / 

When diagonalising these matrices, the fermion mass hierarchies and mixing are well re- 
produced for appropriate values of e, e ~ ~ 0.04. In general, in the models appearing 
in the literature, the relative flavour charges in Table [T] and thus the exact structure of 
the mass matrices are determined by the GUT multiplet structure (and the requirement 
that particles in the same GUT multiplet have the same charge). Nevertheless, in all cases, 
the observed fermion hierarchies require smaller charges for the operators of the higher 
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generations (zero for the top Yukawa mass terms, but also for the bottom and tau in a 
supersymmetric model with large tan/?). This implies that, independently of the specific 
flavour and GUT structure of the theory, and unless extra fields with a non-zero flavour 
charge are involved in the generation of R-violating couplings [22], operators that contain 
fields of the third generation should be naturally larger. 

One has also to worry about the overall suppression of the R-violating couplings with 
respect to the dominant Yukawa ones. However, this may arise either from a small tan (3 in 
supersymmetric models, from the form of the Kahler potential, or from additional, model 
dependent, features of the theory that may involve extra fields and symmetries. 

Along these lines, one may also understand how it could be possible to only have domi- 
nant JJ^DjDk operators. The obvious step, to also ensure the absence of any unacceptable 
proton decay, is to first eliminate lepton-number violating operators by imposing a lepton 
triality, under which the fields transform as 

Z3 : (g, U, D, L, E, H^, H2) (1, 1, 1, a, a^, 1, 1). (3.1) 

This allows only the baryon-number-violating operators and the mass terms, while forbid- 
ding lepton-number-violating one^. In this construction bilinear R-violation would also be 
disallowed. 

To allow only lepton-number violating operators, we could work instead with a baryon 
triality, such as in [9] 

Z3 : (Q, U, D, L, E, Hi, H2) (1, a, a). (3.2) 

Such a baryon triality would allow for bilinear R-violation. However, one may also envis- 
age different structures where the symmetries forbid an LH2 term while allowing trilinear 
lepton-number violating operators. An example of this is given by 

Z3 : (g, f/, D, L, E, Hi, H2) ^ (1, a, 1, 1, 1, 1, a^). (3.3) 

It is interesting to observe that in this case the term fiHiH2 would also be forbidden. 
This is due to charge correlations that arise from the above requirements, plus the need 
to allow Yukawa couplings that generate fermion masses. In this case, the /i-term would 
have to arise either radiatively [26], or through the Kahler potential [27]. The n term 
could also be generated within the framework of the NMSSM [28], via a singlet field with 
appropriate charge; in which case a term SLH2 would also be allowed. Baryon number 
violating operators would be allowed at subdominant orders, due to a term SSUDD which 
is significantly suppressed; moreover, this is not the complete picture, since to explain 
fermion mass hierarchies one would have to introduce flavour dependent charges, which 
could further suppress R-violating operators, particularly for the lighter generations that 
are dangerous for proton decay (see discussion below). 

*A flavour-dependent generalisation of this symmetry has been discussed in [10]. In that scenario, 
consistent solutions were found containing only a subclass of operators violating lepton number (LLE) 
and baryon number (UDD). Thus it is possible to have both lepton and baryon number violation without 
disturbing proton stability. 
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Prom the above, it is clear that whether bihnear or trihnear R- violation dominates is di- 
rectly linked to the symmetries of the underlying theory, and phenomenological information 
would be a valuable probe of this symmetry structure. 

Would these considerations be sufficient to understand the structure of the R-violating 
operators on the basis of positive experimental results? As already discussed, even in 
the case of one dominant operator, for fermions in the basis of current eigenstates, mixing 
effects will induce non-zero coefficients for related operators in the basis of mass eigenstates. 
These will be suppressed by the mixing parameters with respect to the dominant operator, 
but will not be zero, and this may affect phenomenological and cosmological predictions. 
We should also keep in mind that experiments only provide information on the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix \/Ckm _ y^\y^ |^29], and that one can 
construct theoretical models where the left quark mixing is in either the up or the down 
sector, or both. Similarly, lepton mixing comes from the product of matrices of charged 
leptons and neutrinos, with the additional complication that, for the latter, we have the 
possibility of both Dirac and Majorana mass terms (the recent neutrino data indicate the 
existence of neutrino masses and contain the possibility that right-handed neutrinos do 
exist). For instance, in the above mass matrices from [25], the quark mixing is given by 



Due to this mixing, an R-violating operator is in fact a sum of terms. For instance, the 
mixing matrices above would lead to the following interesting mixings: 



These mixings are particularly important since the dominant couplings here have massive 
final states. As we shall see, mixing also opens up for final states forbidden by the gauge 
symmetry of the couplings. However, more generically, we observe the following: 

(i) The right-handed quark mixing (relevant for IJ and /)) is essentially not constrained 
by the data. Therefore, in a model with left-right asymmetric mass matrices, one 
could also imagine a theory with a minimal mixing in the right-handed sector, in 
which case a dominant U^DiDj flavour would be the only relevant one, and a gravitino 
with rriQ < Mt would be essentially stable. 

(ii) For the left quark mixing (relevant for Q), we know the numerical values from Vckm 
(where for instance the 2-3 mixing is a factor of ~ 0.04). Thus, a coupling X'L^Q^D^, 
would in principle also imply the coupling 0.04A'L3(52-D3. 

(iii) The left lepton mixing (relevant for L) is constrained by the lepton data (large 1-2 
and 2-3 mixing, and small 1-3 mixing). 

We see that there are several flavour choices that can lead to significant effects on the decays 
under discussion, particularly in the cases where the available phase space is limited. This 
will be explored in the next Section. 




^U2DiDj + e^UiDiDj + ... 



(3.4) 
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4 Flavour Effects in Gravitino Decays 



4.1 Flavour Effects for LLE Operators 

From the nine R-violating LLE operators, six can potentially give rise to both loop and 
tree-level decays (a common flavour in E and one of the L fields is needed to form the 
loop): 

L2,2,LlEi, Z/i^3L2i?2; L12L3E3, (4-1) 

while three have only three-body decays 

L2L3E1, L1L3E2, LiL2E^. (4.2) 

The cases with a muon or a tau in the loop were discussed in [2]. For an electron in the loop, 
the photonic gravitino decays are very suppressed due to the electron mass, and the tree- 
level decays dominate unless the gravitino becomes extremely light. This is demonstrated 
in Figure [31 where we plot Amax, the maximum allowed coupling, versus the gravitino mass, 
assuming a common slepton mass of 200 GeV. In doing so, we demand that: 

(i) there is one dominant coupling, 

(ii) the gravitinos can be dark matter, with a lifetime of at least 10 times the current age 
of the universe and that 

(iii) photon production from the gravitino decays, as calculated by Eq. (12. ip . is consistent 
with the bounds on the photon spectrum given in Eq. (12.41) . 

We see that while the photon flux from two-body loop decays puts strong bounds on 
couplings that lead to loops with muons (blue) or taus (red), the couplings with electron 
loops (green) are only affected by the three-body decay photons down to very small gravitino 
masses. For the couplings (14.21) with no loop diagrams the bounds are thus correspondingly 
weak, and follow the bound for L2L1E1. As expected, the neutrino flavour has no effect on 
the bounds from the radiative decay, so results for e.g. L1L2E2 and L3L2E2 are virtually 
identical, save for minute differences near the slepton threshold. 

It is also interesting to note that in the terms with only three-body decays in (14. 2p . 
there is always the possibility for tau production in the final state. Indeed, for L1L2E3 an 
SU(2) singlet r is always produced if kinematically allowed, while for L2L3E1 and L1L3E2 
an SU(2) doublet r is produced, unless the gravitino mass becomes lower or comparable 
to the tau. In this case the factor Li^2-t'3 would only contribute to the tree-level decay 
via the z/^e or z/^/i term. This is observed in Figure [3l where in the three-body dominated 
region, bounds on e.g. L1L3E3 are stricter than the bounds on L1L2E2 and L2L1E1, due 
to the extra photons from the tau decay. One can also notice that the bound on L2L1E1 
is slightly better than on L1L2E2] this is due to more bremsstrahlung from electrons than 
from muons in the final state. 
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Figure 3: Maximum value Amax of R-violating couplings versus gravitino mass, for LLE 
operators. Bounds shown as dashed lines are when considering radiative loop decays only, 
solid lines include photons from three-body decays. The sparticle masses are 200 GeV. 

4.2 Flavour effects in LQD operators 

Out of the 27 R-violating LQD operators, only the following nine can potentially give rise 
to both loop and tree-level decays (a common flavour in Q and D is needed to form the 
loop): 

Ll,2,zQlDl, ^1,2,3^2^2, ^1,2,3^3^3, (4.3) 

while the remaining 18 have only three-body decays. 

In Figure H] we show a comparison of the partial lifetime for the loop and tree-level 
decays for the second and third generation. We choose L3, but this has little significance. 
Comparing to the results for the LLE operators in [2], we observe that with a h quark 
instead of a r in the loop, radiative decays dominate over the three-body ones for a sig- 
nificantly wider range of gravitino masses, up to 40 GeV, for the same sparticle masses 
(200 GeV). This arises both due to the higher fermion mass in the loop, but also due to 
the two bottom masses in the final state. The coupling LiQ^D^ gives rise to either iitb or 
Uibb, and the first term is forbidden by phase space up to high gravitino masses, which can 
be seen as a bump in the L3Q3D3 three-body lifetime near thresholdj^ 

As in the previous subsection we can put constraints on the couplings of the LQD 
operators from gamma rays and gravitino lifetime, as a function of the gravitino mass. The 
resulting bounds for operators with both loop and tree-level decays are shown in Figure [5l 
Due to the increased dominance of two-body decays compared to the pure lepton operators, 

^ Close to the bb threshold at ^ 10 GeV hadronization effects will become important for the three-body 
decay. This is not considered here as the two-body decay clearly dominates in this mass range. 
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Figure 4: Comparison of partial lifetime versus gravitino mass for two-body loop decays and 
three-body tree-level decays for the L2Q3D3 and L^Q^D^ couplings. Couplings have all been 
set to X' = 0.001. 



we have even stronger coupling bounds, in particular for the LiQ^D^ couplings, and there is 
now also a significant constraint on the first generation loops, i.e. LiQiDi, for low gravitino 
masses. 

If the mass of the gravitino is close to the lightest possible meson for one particular op- 
erator, we may no longer neglect hadronisation effects from the formation of single mesons, 
as opposed to the QCD jet interpretation of the quarks in the three-body decay. In the 
simplest case we would have a two-body final state with a lepton and a meson, such as 
a pion or a kaon, or even heavier mesons if allowed by the structure of the R-violating 
operator and the mass of the gravitino. For instance, the operator L^QiDi will lead to 
T'K^{Tud) or UT-TT^luT-dd), and similar considerations hold for other fiavour combinations. 

However, since the decay into single mesons is only relevant for low gravitino masses, 
this issue can be neglected for operators allowing loop decays. This is because the constraint 
from the loop decay to photon and neutrino is in all lot more stringent than the 

constraint arising from the decay into mesons at these gravitino masses. 

For operators not permitting loop decays the situation is different, but in the cases with 
light mesons the resulting gamma ray constraints are so weak that other constraints on the 
couplings are more important [12]. Thus the only cases where decays into single mesons 
can have some effect are for operators which give heavy mesons, i.e. B or D mesons. In 
these cases there can be small modifications on the constraints in the range of gravitino 
masses close to the heavy quark masses, but the nature of heavy quarks as kinematically 
equivalent to their corresponding mesons should limit this effect. 
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Figure 5: Maximum value Amax of R-violating couplings versus gravitino mass, for LQD 
operators. The sparticle masses are 200 GeV, except where indicated. 



4.3 Flavour effects in UDD operators 

In this case, we only have tree-level gravitino decays, and of particular interest is the 
possibility of gravitino decays via a dominant U^DjDk operator. For light gravitinos, since 
top production in the final state is kinematically forbidden, decays will arise either 

(i) due to t-c mixing and other possible mixings, 

(ii) or from four-body final states with an off-shell top quark and possibly an off-shell W , 
and with at least one massive final state particle (b-quark). 

The first case is expected to dominate since the second is very suppressed, and the dominant 
decay width should be a function of the right-handed tJ^ — U2 mixing. In this case Amax 
can be large, with interesting phenomenological implications that we discuss in the next 
Section. Another interesting feature of mixing is that it opens up gravitino decay channels 
that were disallowed by the flavour structure of the superpotential, e.g. the possibility of 
two h (or h) in the final state. 

Both of these effects are shown in Figure [6l where we plot the partial lifetime for a 
selection of gravitino decay modes as a function of gravitino mass. We assume a dominant 
coupling \'li2 = 1-0 that for low gravitino masses relies on mixing effects in the decays. 
The coupling is chosen large to minimize the lifetime. We illustrate a possible realization 
of mixing with the mixing matrices in Eq. ([3]), taking e = 0.04 and e = 0.20. As expected, 
it is the t-c mixing that dominates gravitino decays at low masses, and the gravitino is long 
lived enough to be dark matter for a large range of masses. Only for gravitino masses above 
200 GeV, when the top production threshold has been passed with good margin, do the 
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top channels dominate and the gravitino becomes disallowed as a dark matter candidate 
due to its short lifetime. 




Figure 6: Partial lifetime versus gravitino mass for gravitino decays into various quark final 
states with — ^■^> ^ — 0-04 and e = 0.20. All squark masses have been set to 500 GeV. 



We find that changing between the three possible Ag^-^ couplings only changes the relative 
importance of the down type quarks in the gravitino decay, e.g. for Ag^g = 1.0, G cdb 
is the dominant decay channel for low masses. Among the channels that are closed in 
the absence of mixing, we only show decays to b quark pairs. The lighter pairs have very 
similar behaviour to other light quark pairs. We see that the probability of two b quarks in 
the final state is negligible because of the large suppression due to mixing and kinematics 
when compared to other decay channels. Other choices for the mixing matrices only change 
the relative importance of the different decay channels, not the behaviour as a function of 
gravitino mass. 

In Figure [7] we also show the resulting bounds on the A" couplings when considering the 
photon spectrum as in the previous subsections. We notice that the first two generations 
have a log-linear behaviour in terms of the gravitino mass, with equal slopes. The difference 
in scaling is due to different squark masses. With the same squark mass, the two curves 
would be indistinguishable. The importance of mixing effects are again shown for the 
X'^i2 coupling: the opening up of decays through mixing strengthens the bounds on that 
coupling. 

In general, due to the structure of these operators we produce either three quarks or 
three anti-quarks. If there is sufficient phase space, one could imagine that we can end 
up with two-body final states with a baryon and a meson for very light gravitino masses. 
However, we need to keep in mind that the lightest fiavours for A", in particular A'/^2 
and A'/^3, are extremely constrained from double nucleon decay and neutron- antineutron 
oscillations respectively [30]. 
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Figure 7: Maximum value Amax of R-violating couplings versus gravitino mass, for U DD 
operators. The squark masses are as indicated. 

5 Prospects for R-violation in colliders 

For R-violating couplings above 10~^ for 100 GeV sparticle masses, and with a scaling 
that for most operators is a simple proportionality relation, one would expect interesting 
signatures like multi-lepton and/or multi-jet events in the final state of sparticle production 
in a collider. Depending on the flavour of the R-violating operator, the nature of the NLSP, 
and the respective Amax that we found in the previous section, one would generically expect 
either: 

(i) possible observable single superparticle productions, if A can be sufficiently large [5, 
6,31], 

(ii) MSSM production of sparticle pairs followed by R-violating decays of the NLSP, for 
the flavours where Amax is smaller than ~ 10~^, or 

(iii) no R-violatig decays of the NLSP inside detectors for very small Amax (smaller than 
~ 10^^), with some cross-over region where displaced vertices could be observed. 

From the results shown in the previous Section, the observation of single sparticle 
production at the LHC is almost entirely excluded in the gravitino dark matter scenario 
for operators that give loop decays with second or third generation loop-particles, due to 
the strict bounds from gamma rays. For dominant three-body decays the same conclusion 
holds unless the gravitino mass is small {rriQ ^ 10 GeV). Thus the possible astrophysical 
observation of gravitino decays will have important consequences for LHC expectations, 
and vice versa. It is worth noting that this conclusion, for the case of dominant two-body 
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decay, is only weakly dependent on the assumed masses of the other sparticles, as can 
be seen from the insensitivity of Xmax to large changes in intermediate sparticle mass, see 
Figure 5 of [2]. 

For no operator do the constraints considered here eliminate the possibility of seeing 
R-violating decays in colliders, but the LiL^E^i and LiQ^D^ operators allow only a very re- 
stricted coupling range for intermediate to high gravitino masses. Indeed, even for couplings 
of the order of 10^^ it should be possible to detect the R-violating NLSP decays [31]. 

The discovery of supersymmetry at the LHC and the reconstruction of a neutralino 
NLSP has been shown to be possible [32] at the same level or better than for R-parity 
conserving scenarios when one considers the lepton number violating operators. This is 
due to the numerous leptons expected in the final state. However, for the case of a UDD 
operator, assumptions have to be made, either for the production of additional leptons in 
the event from cascade decays, or for heavy flavours that can be tagged. The heavier the 
flavours, the better the detection prospects due to flavour tagging or top reconstruction. 

Decays of the NLSP are highly dependent on the combination of NLSP flavour and 
dominant R-parity violating operator flavour. If these flavours are the same, all NLSP 
decays should be rapid two-body decays if kinematically allowed. In other scenarios, three 
or even four-body decays are the leading decays, with resulting suppression due to phase 
space and heavy virtual particles. For an U^DjDk operator, we have the following particular 
implications: 

(i) we have the possibility of large R-violating coupling with resulting resonant single stop 
production [34] or single gluino top production [35]. Moreover, for a large Ag^^Afg^ 
product, one may observe interesting signatures in single top-bottom production [36]. 

(ii) if the NLSP is a neutralino with a mass larger than the top, it should have a rapid 
three-body decay with a top in the final state, on the other hand, if the neutralino 
is lighter than the top, then it should decay via either subdominant operators or 
mixing effects, which may well enhance its decay rate enough for it to decay within 
the detector, giving a displaced vertex. 

Taken together this would imply the interesting possibility of sparticle production via 
one operator, and decay via a different one. 

6 Conclusions 

We have studied slow gravitino decays originating from lepton or baryon number violating 
operators in R-violating supersymmetry, focusing on the fiavour structure of the theory. 
We found that the dominant decay modes, and thus the final state products are particularly 
sensitive to the hierarchies of R-violating operators and exhibit distinct correlations, which 
we have analysed. Already the dominance of trilinear R-violating couplings over bilinear 
modes implies the presence of symmetries that, among others, have interesting implications 
for the /i-term. 
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A more detailed study of the flavour dependence of the operators has determined the 
ratio between (i) the tree-level gravitino decays to three fermions and (ii) the two-body loop 
decays into a photon and a neutrino, which in turn puts strong bounds on the maximal value 
of the allowed R-violating couplings. Bounds from photon spectra are much stricter than 
the ones from the requirement on the gravitino lifetime, and thus strongly constrain the 
respective operators, particularly LiL^E^ and LiQ^D^ that involve a r and a bottom-quark 
in the loop. On the other hand, for operators without photonic decays larger coupling 
constants are possible, particularly in the case of phase space suppressions due to the 
presence of heavy fermions in the final state. Moreover, mixing effects turn out to be 
crucial in correctly deriving the relative magnitude of the various contributions, and play a 
significant role for decay modes with phase space suppression and particularly for the ones 
generated by U^DjDk. 

In all cases, the bounds on the R-violating couplings from the cosmological requirements 
are compatible with visible signatures at colliders, which can vary from single superparticle 
production (for flavours where a larger coupling constant is allowed) to MSSM production 
and R-violating decays (for the smaller couplings). Particularly for the operator flavours 
that would lead to predominantly photonic gravitino decays, giving strong constraints on 
the couplings, interesting event properties such as vertex displacement might be expected. 



A Photonic Gravitino Decays 

The photonic decays of the gravitinos have been calculated in [2] . The rate for the radiative 
decay G — >^ 7Z/ with the loop fermion / is given by 

20487r4 M2'-^' ' ^ ' 

where Mp = (SttGat)"^/^ = 2.4 x 10^^ GeV is the reduced Planck mass, anc^l 

1 , ,0 2, ,0 1, 



= -\c,\' + -\C2\' + ^Reiclc), (A.2) 
with 

ci = 2[(m| - m^- + m})d,^^ + ASj'^], 

C2 = 2 [mjci"^) + mj- Cj") + A^f ] , (A.3) 

where, in the notation of LoopTools [37,38], we have 

C^"^ = Co(m|, 0, 0, mj, mj, mj), 
Cq^ = Co(m|, 0, 0, m^, mj, mp. 



^Here, we correct a minor error in that paper due to a misprint in the gravitino spin-sum taken from [33], 
where the sign in Eq. (4.31) should be {-^ + 7713/2)- 
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AeJ^^ = 2Eo(m|,mJ-,mJ) -5o(0,mJ-,mJ) -Eo(0,m5,mJ), 

A5f = 5o(m|, mj-, mj) - Bo(0, mj-, mj). (A.4) 

The Co are three-point functions corresponding to Fig. [T](a) and (b), whereas the ASq are 
finite differences of two-point functions. 



B Three-body Gravitino Decays 

The three-body decays of gravitinos have been calculated in [4], where extensive analytic 
formulas were derived. Here, we only comment on the spin summed squared amplitudes, 
and refer to the original paper for the full computation. 

The full squared amplitude (summed over spins) for the gravitino decay G Ujk^^dj^k 
is the sum of three individual squared amplitudes plus three interference terms. These arise 
since the gravitino can couple to all the particles involved in the R- violating operator. Then, 
for the case where the gravitino couples to a neutrino and a sneutrino, one has 



Hjk 



X 



m-jk j , 



(B.l) 



where m^^ = {pj + PkY, with pj^k the four-momenta of the respective particles. The 
remaining squared amplitudes Mb^c, where the gravitino couples to the charged lepton of 
the doublet and the singlet charged lepton respectively, are given by the same formula, when 
substituting the appropriate flavours in the vertices and the propagator. The interference 
terms are of the form 



2Re {MaMl 



1 



''ijk 



3 M2(m|^ - mlj{mj^ -mj ) 
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(B.2) 



For LQD operators the results are similar, and found by replacing the SU(2) doublet 
fleld L by Q, and the SU(2) singlet E by Z), and summing over colours. For tJDD operators 
we also have similar amplitudes and interference terms. Again the contributions can be 
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read off from Eq. (IB.ll) and (]B.2I) . modulo colour and symmetry factors that arise from the 
possibility of two identical particles in the final state. 
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